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Abstract. The existence of magnetohydrodynamic mean-field a 2 -dynamos 
with spherically symmetric, isotropic helical turbulence function a is related 
to a non-self-adjoint spectral problem for a coupled system of two singular 
second order ordinary differential equations. We establish global estimates for 
the eigenvalues of this system in terms of the turbulence function a and its 
derivative a' . They allow us to formulate an anti-dynamo theorem and a non- 
oscillation theorem. The conditions of these theorems, which again involve a 
and must be violated in order to reach supercritical or oscillatory regimes. 



1. Introduction 

One of the simplest models of a magnetohydrodynamic dynamo is a mean-field 
a 2 -dynamo with spherically symmetric, isotropic turbulence function (a-profilc) 
a(r). Models of this type were among the first that were capable to explain the 
turbulence-based inverse energy cascade leading to non-decaying dynamo regimes 
with self-sustaining magnetic fields (see |KR80] ). Recently, a 2 -dynamo models were 
used for detailed studies of polarity reversal processes of magnetic fields (see, e.g., 
SG05 , SGGX06j), as suggested by paleomagnetic data of the Earth's magnetic 
field. Deep insight into the reversal dynamics could be gained by analyzing the 
highly non-linear back-reaction-based induction processes by a time foliation method 
built over a series of instantaneously linearized (kinematic) auxiliary setups. For 
the latter, the a-profiles were assumed to be fixed radial functions so that the well- 
known decomposition of the magnetic field into poloidal and toroidal components 
with subsequent expansion in spherical harmonics could be used. As a result, one 
arrives at a set of eigenvalue problems, indexed by the degree i€N = {l,2,...}of 
the spherical harmonics, for pairs of coupled linear ordinary differential equations 

in £2(0, 1) © £2(0, 1) subject to the boundary condition 
(1.2) ( {9r + l)yi )(l) = 



V2 

(see [KR80| Section 14.2], [SG03| (6)-(8)]). In the following, for brevity, we call the 
eigenvalue problem (jl.ljl . (|1.2|) dynamo problem. 
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In view of the time separation ansatz B(x, t) = exp(Ai)b(x) of the corresponding 
magnetic field modes, one naturally distinguishes between decaying or sub-critical 
modes (Re A < 0) and amplifying or supercritical modes (Re A > 0), as well as 
between oscillatory modes (ImA ^ 0) and non-oscillatory modes (ImA = 0). The 
physically relevant self-sustaining dynamo configurations are mainly defined by a 
few supercritical modes, whereas possible polarity reversals of the magnetic fields 
are closely related to the existence of oscillatory modes close to criticality (Re A m 0, 
ImA 7^ 0) (see |SG05j . [SGGX06] ). For a deeper understanding of the physical 
interplay between the plasma- (or conducting fluid-) based a-profile and the dynamo 
dynamics, knowledge of the spectral properties of the dynamo problem JO}, (H3) 
is of utmost interest. 

In contrast to the large number of general results for kinematic fast dynamo 
problems (see, e.g., [AK98[ Chapter V] as well as CLMS95]) and despite the sim- 
plicity of the physical model, no proven analytic information for the spectrum of 
the kinematic mean-field dynamo problem (jl.ip . (|1.2[) seems to be available, except 
for the case of constant a where the spectrum is real. For non-constant functions 
a, only numerical calculations for eigenvalues were performed (see, e.g., |SG03| ). 
However, for a non-symmetric spectral problem like (| 1 .2[) . numerical compu- 

tations are prone to be unreliable; a convincing example for this is a non-normal 
7x7 matrix due to S.K. Godunov for which numerical algorithms yield results far 
away from the true eig envalues (see, e.g., [God98| . |HP051 Example 5.2.5]). 

Other attempts to attack the above dynamo problem include replacing the phys- 
ical boundary conditions by the so-called idealized boundary condition 

(1.3) © (1)=0 ' 

assuming that the helical turbulence function a is strictly positive, and/or analyzing 
the simpler monopole modes, that is, the case I = 0, (see, e.g., [GS03] . [GK06 , 
KGS09 ); all these restrictions are not desirable from the physical point of view. 

The aim of this paper is to establish analytic enclosures for the spectrum of the 
dynamo problem (11.11) . (jl.2l) in terms of the function a, not making any of the 
above simplifying assumptions. We show, in particular, that the eigenvalues lie in 
a horizontal strip around the real axis and we establish an upper bound for their 
real parts. More precisely, every eigenvalue A of (jl.ip . (|1.2j) satisfies 

(1.4) |ImA| < ||a'||, Re\<p e 

where the constant pg depends on ||a|| and ||a'||, as well as on the smallest eigen- 
values Ai (Z) and Ai(oo) of the Bessel operator —d 2 + 1(1 + l)/r 2 with boundary 
condition y'(l) + ly{l) = and y(l) = 0, respectively (compare (11.21) ). Here, for a 
continuous function q on [0, 1], we denote by \\q\\ := max rg [ .i] \q(r)\ the maximum 
norm of q. The estimate for the real part in (|1.4I) yields, in particular, a so-called 
anti-dynamo theorem: if 

II II? Il a llll a 'll x / X 

IH| 2 + 7= < Ai(oo), 

V A H*) 

then the dynamo operator has no eigenvalues in the (supercritical) right half-plane 
(compare the more geometrical anti-dynamo theorems for kinematic fast dynamos 
listed, e.g., in [AK981 Chapter V, §3]). In addition, we establish conditions ensuring 
that a particular eigenvalue remains on the real axis when a ^ 0, which could be 
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called a local non-oscillation theorem. Our results also cover other boundary condi- 
tions such as the idealized case (|1.3|) . and they do not require a to be sign definite. 
The methods we employ stem from the perturbation theory of linear operators and 
from the spectral theory of block operator matrices (see, e.g., |Kat76| . |Tre08] ) . 

In the following we give a brief outline of the paper. In Section[2]we introduce a 
family of linear operators Ag in £2(0, 1)©£2(0, 1) so that we can write the boundary 
eigenvalue problems (|1.2j) and (jl.lj) . (|1.3|) as spectral problems (Ae — A)y = 

with 6 — I and 6 — 00, respectively. The operators Ae are block operator matrices 
of the form 



Here Ag and Ag a are linear operators in L 2 (0, 1) given by the Bessel (type) differ- 
ential expressions — d 2 + l(l + l)/r 2 and — d r ad r + al(l + l)/r 2 , respectively, and the 
boundary condition x' '(1) + 6x(l) = for 9 G [0, 00] (note that 8 = 00 corresponds 
to cc(l) = 0). In Section[3]we investigate the entries of the block operator matrices 
Ag in detail and collect all the properties needed in the subsequent sections, e.g., 
to show that Ae defines a closed linear operator and that its spectrum consists only 
of eigenvalues of finite algebraic multiplicities with no finite accumulation point. 

Section [¥] contains our first main result, the eigenvalue enclosure in Theorem l4.6l 
This first estimate is based on decomposing Ae into a lower triangular block oper- 
ator matrix plus a bounded part containing only the right upper corner a: 



This decomposition, together with a Neumann series argument, enables us to iden- 
tify a region in the complex plane in which the eigenvalues must lie. In addition, 
we establish bounds on ||a|| and ||a'|| such that an eigenvalue of one of the diag- 
onal elements of Ae, that is, an eigenvalue when a = 0, remains real for a ^ 
(see Propositions 14. 101 and 15. 12|) . Physically speaking, this means that a particular 
mode remains non-oscillating. 

Section [5] contains our second main result, the eigenvalue enclosure in Theo- 
rem ^. 101 Here we use a quasi-similarity transformation of Ae (with the unbounded 

1 /2 

operator We = diag (A e ' , /)) such that the transformed operator Be is a bounded 
perturbation of a self-adjoint operator: 



This allows us to conclude that the eigenvalues of Ae lie in discs of radius ||a'|| (the 
norm of the perturbation Te) around the eigenvalues of the self- adjoint operator Se] 
in particular, the imaginary parts of the eigenvalues of Ae are bounded by ||a'|| and 
their real parts are bounded by max a(Se) + \\ct'\\. 

In Section [5] we compare the two different eigenvalue enclosures of Sections 0] 
and[5j It turns out that, generically, a combination of both estimates yields the best 
result: it consists of the uniform estimate for the imaginary parts of the eigenvalues 
from Section [5] and of the upper bound for the real parts of the eigenvalues from 
Section [4] (see (|1.4j) and Figure |4j . At the end of Section |6] we summarize the 
physically most relevant result in the form of an anti-dynamo theorem (see (|6.3j) ). 
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In Section [7] we illustrate our results by some examples. They include the case 
of constant a where the eigenvalues of the physical dynamo problem are only given 
implicitly as solutions of an equation involving four different Bessel functions. We 
also consider the particular non-constant function a for which dipole-dominated 
oscillatory criticality was first found numerically by F. Stefani and G. Gerbeth (see 
SG03]); in their computations they obtained special a-profiles such that the right- 
most eigenvalues pass from Re A < to Re A > with Im A 7^ first for the dipole 
modes (1 = 1) and only afterwards for quadrupole and higher-degree modes (I > 1). 

The following notation is used throughout the paper. In the Hilbert space £2(0, 1) 
the scalar product and norm are denoted by (•, •) and || ■ || , respectively. By ACi oc (-0 
we denote the space of locally absolutely continuous functions on a (sub-) interval 
/ C [0, 1] (that is, the space of functions that are differentiable (Lebesgue-) almost 
everywhere on / and have locally Lebesgue-integrable derivative on /). For a lin- 
ear operator T we denote by T>{T) its domain, by cr(T), a p (T) 7 p(T) its spectrum, 
the set of its eigenvalues, and its resolvent set, respectively; if T is bounded, we 
denote by ||T|| its operator norm. For a continuous function q G C([0,1]), we use 
the same symbol q to denote the (bounded) multiplication operator by the func- 
tion q in £2(0, 1); in this case, the operator norm is given by the maximum norm 
\\q\\ = maxy-gro^i |g(r)|. For details on linear operators in Hilbert spaces we refer 
the reader to |Kat76l Chapter III, § 2, 5, 6], [RSlOl Chapters VI, VIII], and [RS78] , 

Finally, we would like to thank Frank Stefani and Gunter Gerbeth (Research 
Center Dresden-Rossendorf) for useful discussions on the physics of dynamos and 
for sharing their numerical data with us. Special thanks also go to Markus Wa- 
genhofer (now Psylock, Regensburg) who performed various numerical calculations 
while this manuscript was being prepared. 

2. Operator model 

In this section we introduce a family of linear operators Aq in the product Hilbert 
space £2(0, 1) © £2(0, 1) which describe the system of differential equations (II. ip 
with physical boundary conditions (11.21) and idealized boundary conditions (11.31) . 

To this end, we first associate linear operators in the Hilbert space £2(0, 1) with 
the two differential expressions 

(2.1) r := -d; + ^ , r Q := ~d r a(r)d r + a(r) — — — , d r := — , 

occurring in (jl.ip . Here and in the following, we always suppose that I G N = 
{1, 2, . . . } is fixed and that a: [0, 1] — > M is a continuously differentiable real-valued 
function, a € C 1 ([0, 1]). 

The classical Bessel differential expression r appears twice on the diagonal in 
(ll.ip . but it is subject to two different boundary conditions at r = 1 in (|1.2p . 
Therefore we define a family of operators Ag, 6 € [0, 00], by means of the boundary 
condition y'(l) + 9y(l) = 0, so that 9 = I yields the first order boundary condition 
for yi and 9 — 00 yields the Dirichlet boundary condition for y% in p.2[) . 

Definition 2.1. Let 9 £ [0, 00]. Define the linear operator Ag in £2(0, 1) by 
V(A e ) := {x g £ 2 (0, 1) :x,x'e AG oc ((0, 1]), tx g L 2 {0, 1), x'(l) + 9 x{\) = 0}, 

(Agx)(r) := (rx)(r) = —x"(r) + - ~^ — -x(r) for almost every r G (0, 1]. 
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Remark 2.2. Every x £ T>(Ag) automatically satisfies linv^o x(r) — 0. This follows 
easily if one calculates the inverse Ag 1 and checks that lim r \ i o(A g ~ 1 f)(r) = for 
/ G L 2 (0, 1) (see, e.g., |AHM07l Lemma 2.6 and its proof]); note that a fundamental 
system of ry = is given by the simple functions x\{r) — r~\ X2(r) = r +1 , r£ (0, 1]. 



The differential expression r a in (|2.ip is not classical since its leading coefficient 
a may change sign and may hence also have singularities in the interior of the 
interval [0, 1]. Since r Q is only located in an off-diagonal corner in (jl.l[> . the following 
definition is sufficient for us. 

Definition 2.3. Let 9 £ [0,oo]. Introduce the linear operator Ag tCi in £2(0, 1) by 
V(A e , a ) :=V{A e ), 

(Ag a x)(r) :— (r a x)(r) = — (ax')'(r) + a(r) — \ ^ x(r) for almost every r £ (0, ll. 

r z 

Remark 2.4. In Proposition l3.6l below we show that Ag a x £ L 2 {Q 1 1) for x £ T>(Ag), 
that is, Ag a is well-defined. 

Now we are able to formulate the dynamo problems (|1.1[) . (|1.2p and (|l.ip . (|1.3p 
as spectral problems for a linear operator acting in the product Hilbert space 
£2(0,1) © L 2 (0, 1). Here the case 9 = I corresponds to the physical boundary 
condition (II. 2[) , while 9 — 00 corresponds to the idealized boundary condition (| 1 . 3[) . 

Proposition 2.5. Lei £ [0, 00]. Define a linear operator Ag in £2(0, 1)0^2(0, 1) 
by the block operator matrix 

(2.2) Ag:=( ~ A . B a ), V{Ag):^V{A e )®V{A OQ ). 

Then the boundary eigenvalue problems (jl.ip . (jl.2p and (jl.ip . (|1.3p can 6e written 
equivalently as 

(Ag-X)y = Q, y£V(Ag), 
for 9 = I and 9 = 00, respectively. 

Proof. The claim is evident from the definitions of the operator Ag and of its entries 
Ag and Ag <a for 9 = I and 9 = 00. □ 

3. Auxiliary results 

In this section we study the properties of the entries Ag and Ag. a of the block op- 
erator matrix Ag introduced in the previous section. They are used in the next two 
sections to derive estimates for the eigenvalues of the dynamo problem (jl.ip . (|1.2I) . 

Proposition 3.1. Let 9 £ [0, 00] . The linear operator Ag in Definition ® .11 has the 
following properties: 

i) Ag is self-adjoint and positive; the domain of its square root is given by 

(3.1) P(Aj /2 ) = € L 2 (0, 1) : 16 AC loc ((0, 1]), x', ^gL 2 (0, 1)}, 0g[O,oo), 

(3.2) V(A^ 2 ) ={x G L 2 (0, 1): x£ AC loc ((0, 1]), x' , ^GL 2 (0, 1), *(1) = o}. 

ii) Ag has compact resolvent, and the spectrum o~(Ag) consists of a sequence of 
simple eigenvalues < \i(9) < Xzif)) < ■ ■ ■ tending to 00. 
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iii) Xk(0) is strictly increasing in 8 for every k — 1, 2, ... , and for different val- 
ues 9\, 82 the sequences (Xk{di)) k (A J t(^2)) fe „ 1 interlace; in particular, 

(3.3) < Ai(0) < Ai(oo) < A 2 (0) < A 2 (oo) < ••• , 0e[O,oo). 

iv) For < 81 < 62 < 00 and A > 0, we have 

(Ae 2 +A)- 1 < (A 6l +X)-\ 

Proof. All claims of the proposition are well-known; for convenience of the reader, 
we repeat the main arguments. 

i) The Bcsscl differential expression t (with I > 1) is in limit point case at 
the singular end-point (see, e.g., [AG93, Appendix II, Section 9. IV]); in fact, 
the £2(0, l)-solutions of the differential equation (t — X)x — are spanned by the 
Riccati-Bessel function (see |AS64[ 10.3.1]) 

(3.4) f t (r, A) := rVXj t {rVX) = ^\fr\7^ J l+1/2 (rVx), r £ [0, 1]. 

This implies that the operator Ag as denned above is self-adjoint (see |Wei03[ 
Satz 13.21a)]); moreover, the set 

(3.5) T> (Ag) := {x £ V(Ag) : suppx compact in (0, 1]} 

forms a core of Ag, that is, the closure of the restriction Ag\x> (A e ) coincides with Ag 
(see |Nai68| Section 17.4]). 

The positivity of Ag follows from the identity 

(3.6) ag[x] := (Agx,x) =: (x,x)g, x £ V(Ag), 
where 

(3.7) (x,y) e :=c e x{lj^)+ [ x\r)J[r)dr + 1(1 + 1) / ^IfM dr 

Jo Jo T 

with eg = 8 for 8 £ [0, 00) and Coo = 0. The inner product (x, y)g is defined for all 
x,y £ £2(0, 1) such that x' , y' , x/r, y/r £ £2(0, 1); in particular, it is defined on 
the sets on the right-hand sides of (13. 1[) and of (|3.2I) . 

The domain T>(A g ' ) is the domain V(ag) of the form closure ag of the qua- 
dratic form a e given by (J3U) (see |Kat76[ Theorems VI.2.23, VI.2.1, and Corol- 

1 II 

lary VI. 2. 2]). Hence for every x £ D^Ag ) there exists a sequence (xn)^ C V(Ag) 
such that x n — > x in L 2 (0, 1), n — > 00, and 

(3.8) (Ag(x m x n ), x m x n ^ — (x m x n , x m x n )g y 0, ?n, n > 00. 

The relations (pT5)) . (|X71) yield that a^, - x' n 0, z m /r - x n /r ^ in L 2 (0, 1) 
and, if eg ^ 0, also x m (l) — x„(l) — > for m, n 00. Since convergence in 
£2(0, 1) implies convergence almost everywhere, we can choose ro £ (0, 1] such that 
Xni^o) ~~ ^ x(ro), n — > 00; if eg ^ 0, we can always choose ro — 1. Together with 

x n (r) - x n (ro) = / x'(t)dt, re [0,1], 

it readily follows that the sequence (£n)o° converges uniformly in [0, 1] to x, that x 
is absolutely continuous with x' £ £2(0, 1), and that x' n — > x' in L 2 (0, 1), n — > 00. 
The relation (|3.8p also shows that the functions x n /r form a Cauchy sequence in 
£2(0, 1) and hence also x/r £ L 2 (0, 1). If 8 = 00, then x n £ V(A 00 ) implies that 



SPECTRUM OF MAGNETOHYDRODYNAMIC MEAN-FIELD DYNAMO OPERATORS 7 



x n (l) = 0, n G N, and due to the uniform convergence of (x n ) it follows that 
x(l) = 0. Thus we have proved the inclusions "c" in (|3.ip and (|3.2p . 

In order to prove the converse inclusions "D" in (|3.1I) . (|3.2[) . we equip the set 
He on the right hand side of (13. ip or (|3.2p . respectively, with the inner product 
(•, -)e defined in (|3 . T|) . Then He becomes a Hilbert space and, by its definition, the 
subset V(ag) is a closed subspace of He- Now assume that y £ He is orthogonal to 
T>(ag). Then, in particular, for every x £ T>q(A$) C T>(Ae) = V{ae) with T> (Ae) 
given by (|3. 51) , 



(3.9) = (.T,y )e = (Aea;, y )- 

Since I>o(-^e) is a core of Ag, we conclude that yo £ D{Aq) — T>(Ae) and Aeyo = 0. 
The positivity of Ag now implies yo = 0. This proves that DfAg ) = 2?(ae) = "He- 
ii) In |AG93[ Appendix II, Section 9. IV] it was proved that the spectrum of Aoo, 
and hence of every other operator Ae (see [AG931 Appendix II, Theorem 6.2] or 
[Wei03|, Abschnitt 14.2]), is discrete or, equivalently, Ae has compact resolvent; 
moreover, all eigenvalues are simple. In fact, if we choose a fundamental system 
{<p(-, A),^>(-, A)} of (r — X)x — such that 

¥>(1,A)=0, V(1,A)=1, 
^'(1,A)=-1, V'(1 J A)=0 > 

where ' denotes the derivative with respect to the first variable r, then the corre- 
sponding Weyl-Titchmarsh function m,oo is given by (see |Tit62[ (4.8.2)]) 

The Weyl-Titchmarsh function me for € [0, oo) can be written as 

1 + 6> (A) . _ 

p - m-oo(A) 

(see |Wei03[ Lemma 14.10]). Since the eigenvalues Afc(0) of Ae are the poles of the 
Weyl-Titchmarsh function me, Afe(co) is the square of the fc-th non-zero zero of the 
Bessel function J;+i/2, while Xk(fi) is the fc-th non-zero zero of the function moo — 
for 9 £ [0, oo). For all 9 £ [0, oo], the eigenspace of Ae at an eigenvalue Afc(6>) is 
spanned by the function 



x k .e{r) = fi(r^/\~US)) = J^yfr V^) J l+1/2 (ry/xWj) , 



6(0,1]. 



iii) The function moo is a Nevanlinna function, that is, it is analytic in the upper 
half plane C + , maps C + into C + U R, and is complex symmetric with respect to 
R (see, e.g., [AG931 Appendix II, Section 9.4 and Section VI.59, Theorem 2]); in 
particular, moo is real-valued on R. Moreover, moo is meromorphic, its singularities 
are the simple poles Afe(oo), and it is strictly increasing between two neighbouring 
poles. This and the property that Xk(9) is the fc-th zero of moo — for 9 £ [0,oo) 
imply both claims of iii). 

iv) The claimed inequality is an immediate consequence of Krem's resolvent 
formula (see |Kre46j . |AG93[ Section VIII. 106], or |GM09i Section 9]). □ 



^Note that, in the first edition of ITit62l from 1946, the term —1/2 was missing. 
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Remark 3.2. By Proposition 13.11 i). for 9 — oo the boundary condition x(l) = 

1/2 

carries over from T>(A 00 ) to T){AJ 2 ); it is called an essential boundary condition 
(comp. |Kre47t §7]). For 9 £ [0, oo) the boundary condition x'(l) + 9 x(l) = 

1 /2 

in T>(Ag) disappears in T>(A g ); it is called a non-essential or natural boundary 
condition (see |BS87| Section 10.5, p. 234])). As a consequence, 2?(Ag /2 ) does not 

1 /2 

depend on for 9 € [0, oo), whereas the action of A g does. 

In order to study the problems (jl.ip . (|1.2[) and (jl.l[) . (| 1 . 3[) . we shall need, in 
particular, the eigenvalues of the operators A^ and Af. 

Lemma 3.3. The eigenvalues Xk(9), k — 1,2,... for 9 = I and 9 = oo are as 
follows: 

i) Afc(oo) is the k-th non-zero zero of the function A i— > J;+i/2(v / A) ; 

ii) Afc(Z) is the k-th non-zero zero of the function A M- Ji-i/2(V~X), 
for k= 1,2,.... 

Proof. For = oo, the boundary condition xi(l) = implies that A € C is an 
eigenvalue of if and only if </;+i/2('\/A) = 0. For 9 = 1, the boundary condition 
x',(l) + lxi(l) = implies that A € C is an eigenvalue of A\ if and only if 

(3.10) = ^(Vr J m/2 (rVA)) +1^ J l+1/i {ry/\)\^ = J ; _ 1/2 (VA); 

here we have used the differentiation formula (z l -zji(z)) = z' +1 j;_i(z) for the 
spherical Bessel functions ji (see |AS64[ 10.1.23]). □ 

Remark 3.4. Note that the index 9 in Ag and its eigenvalues Xk(9) only refers 
to the boundary condition x'(l) + 9x(l) = 0. Since the differential expression 
r = —d% + 1(1 + l)/r 2 defining Ag always depends on I, so do all eigenvalues Xk(9) 
(compare Lemma T3.3p . As I £ N is assumed to be fixed throughout the paper, we 
do not indicate this dependence in general. 



For the differential expression r a in (|2.1I) . we use the relation 

r a = a t — a'd r , 

thus relating the operator Ag^ a induced by r Q to the operator Ag induced by r. 
Lemma 3.5. Let 9 6 [0, oo] and let D be the operator of differentiation in £2(0, 1), 

V{D) := 1^2(0,1), Dx:=x'. 
Then the operators DA g and DAg 1 are defined on £2(0, 1) and bounded with 

(3.11) \\DA g 1/2 \\ < 1, ll^- 1 ]]^ ' 



Proof. Proposition 13. H i) and its proof yield that V(A 1 J 2 ) = T>(a$) C T)(D) and 

\\Al /2 x\\ 2 = a5[x] > / \x'(r)\ 2 dr = \\Dxf, x £V{A l J 2 ), 
Jo 

that is, 1 1 DAg 1/2 1| < 1. The claims for DA g 1 now follow from the identity DA g = 
DA g 1/2 A g 1/2 and from the estimate ||A^ 1/2 || < l/y/Xi(9). □ 
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Proposition 3.6. Let 8 G [0, oo]. The linear operator Ag. a in Definition 12.31 is 
densely defined, symmetric and hence closable, and it satisfies 

(3.12) A e . a = aA e - a'D. 

Proof. First we have to show that Ag iCl is well-defined, that is, ax' G ACi oc ((0, 1]) 
and T a x G £2(0, 1) for x £ T>(Ag). To this end, let x € T>{Ag). By definition, this 
implies x' G ACi oc ((0,l]) and tx £ L 2 (0, 1). Since a G C^QO, 1]), it follows that 
ax' € ACi oc ((0, 1]). Further, Lemma [3.51 shows that T>(Ag) C T>(D) and hence 
x' G £2(0, 1). Together with a, a' G C([0, 1]), we obtain 

r a x — arx — ax G £2(0, 1), 

and also the operator identity (|3.12j) . For the symmetry of Ag tCC , it suffices to show 
that (Ag t0l x,x) G M for all x € 2?(Ag iQ ) = T>(Ag). Using the boundary condition at 
r = 1, we see that 

(Ag :Ct x, x) = a(0) lim (x'(e) x(e)) + c e a(l) \x{l)\ 2 



1 



+ / a(r)|x'(r)| dr + a(r) 



1(1 + 1) 



x(r)\ 2 dr, 

Jo '" 

where eg = 8 for 8 G [0, 00) and = 0. Since x € 2?(-Ae) and since the differential 
expression r defining is in limit point case at 0, the limit in the first term on 
the right hand side is real (see jWei031 Satz 13.19]). Because a is real-valued, it 
follows that (Ag, a x, x) G K. □ 

We close this section with some resolvent estimates for the self-adjoint opera- 
tors Ag, which will be used in the next sections. 



Lemma 3.7. For 8 G [0, 00] and X 

1 



\(Ag + A) -1 1 



< 



|Ai(0) + A| 
1 



ImA 



\K /2 (A g 



+ A)- 1 ! 



\\i(8) + X\ 



\Ag(Ag + A) -1 ! 



< 



< 



|A| + A| 



— 00, — Ai(6*)], we have the norm estimates 
if Re A > -Ai(0), (i) 

if Re A < — Ai(0), (ii) 

if |A|<Ai(fl), (a) 

|A| >Ai(0), (b) 
«/ Re A > 0, (1) 



|A| 



|ImA| 
Ai(0) 



if Re A < 0, |A+ Ai(0)/2| > X 1 (8)/2, (2) 

|Ai(&) + A| |A + A lW /2|<A lW /2; (3) 

m eac/i 0/ i/ie i/iree cases, the bounds on the right-hand side define a continuous 
function of X on C \ (—00, — Ai(0)l . 

Proof. The resolvent of the self-adjoint operator Ag satisfies the well-known relation 
(see, e.g., [Kat76| Section V.3.5]) 

l*| s 



(3.13) 



A s e (A e + xy 



sup + 



G[0,1], A€V(A e ). 
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For s = 0, the right hand side of (|3.13|) equals 1/dist (— A, cr(Ag)) and, together 
with the inclusion a(Ag) C [Ai(6*), oo), the equality (i) and the estimate (ii) follow. 

For s = 1/2, it is not difficult to see that the function f(t) :— \ftj\t + A|, 
t £ [0, oo), has a local maximum at t\ = |A|. Hence, if |A| < \i(0), then the function 
/ restricted to a(Ag) C [Ai(#),oo) attains its maximum at t = Xi(9) £ cr(Ag) and 
the equality (a) follows; if |A| > \\{9), then / restricted to a(Ag) C [Ai(#),oo) is 
estimated by its maximum on [Ai((9),oo) attained at t — |A| and (b) follows. 

For s — 1, a short calculation yields that the function g(t) := t/\t + A|, t £ [0, oo), 
is monotonically increasing if Re A > and has a local maximum at t\ — | A| 2 /|Re A| 
if ReA < 0. Thus, if ReA > 0, then the function g restricted to a(Ag) C [Ai(#),oo) 
attains its maximum at oo and (1) follows. If ReA < 0, we note that 

<a = ^<A 1 (0) |A| 2 + ReA • \i(6) < ^ |A + Ax(<?)/2| < 

and therefore in this case the function g restricted to a(Ag) c [Ai(6*),oo) attains 
its maximum at t = \\{9) £ a(Ag) and the equality (3) follows; if |A + Xi(9)/2\ > 
Xi(9)/2, then g restricted to a(Ag) C [Ai(#),oo) is estimated by its maximum on 
[Ai(#),oo) attained at t\ = |A| 2 /|ReA|. Now another short calculation shows that 
g(t\) coincides with the upper bound in (2). □ 

4. The spectrum of Ae and a first eigenvalue estimate 

In the following we present a first perturbational approach to study the spectral 
properties of the operator Ae associated with the dynamo problem. To this end, 
we regard the block operator matrix Ag as a bounded perturbation of its lower 
triangular part. We use this decomposition to show that Ag is a closed operator 
with compact resolvent and to establish estimates for its eigenvalues. 

Theorem 4.1. Let 9 £ [0, oo] and define 

(4.1) Qg:=( 'j 9 J ), V(Qg)=V(Ag)(SV(A QO ), 

(4.2) ^ : =(o o)' 2>(rc)=L 2 (0,l)©£a(0,l). 
Then the block operator matrix 

(4.3) Ag = Qg+K 

is closed and has compact resolvent; its spectrum is symmetric to R and consists 
of isolated eigenvalues of finite algebraic multiplicities with no finite accumula- 
tion point. 

Proof. Since a £ C([0, 1]), the operator 1Z is everywhere defined and bounded. 
Moreover, T>{Qg) = V(Ag) and hence the operator identity (|4.3p holds. 
Next we show that 

(4.4) p{Q 6 ) = p(-Ag) n pi-Ac) £ 0, 

which implies the closedness of Qg, and that (Qg — A) -1 is compact for all A £ p(Qg). 
By (pTP2]) . we have 

Ag, a {A e + A)- 1 = (a(A 8 + A) - aX - a'D^j (Ag + A)- 1 
= a- (aX + a'D)(A g + X)- 1 . 



SPECTRUM OF MAGNETOHYDRODYNAMIC MEAN-FIELD DYNAMO OPERATORS 11 



Since a, a' G C([0, 1]) and D(Ag + A) -1 is bounded for A G p(— Ag) by Lemma l3?5| 
the operator Ag iCC (Ag + A) -1 is bounded for every A G p(— Aq). Therefore, for 
A G p{— Ag) n p(— Aoc), the inverse 



-(Ag + X)- 1 
-(A^ + X^Ag^Ag + X)- 1 -(A^ + XY 



(4.6) (Q fl -A) 

is bounded and everywhere defined. This proves the inclusion "d" in (|4.4[) . Vice 
versa, assume that Xea(-Ag) U er(— Aoo) — (T p (— Ag) U <7 P (— -Aqo). If Ag<7 p (— Aqo), 
then there is an element y 2 G ker(-^4 oc - A) \ {0}. Then (0, y 2 f G ker(Q - A) \ {0} 
and hence A G o- p (Qg) C o-(Qe). If 9 G [0,oo) and A G cr p (— ^4g), then there is 
an element y\ G ker(— Ag — X) \ {0}. Since 9 7^ 00, we have A ^ cr(— A^) by 
Proposition 13.11 iii) and so (yi, (A^ + X)~ 1 Ag :a yi) t G ker(Qe — A) \ {0}, that is, 
A G <7 p (Qg) C cr(Q e ). This proves (j4~4l) . 

By Proposition 13. II ii) the inverses (Ag + A) -1 are compact and we have shown 
that Ag, a (Ag + A) -1 is bounded for all A G p(—Ag) and 9 G [0,oo]. Thus, by 
equation (|4.6[) . the inverses (Qg — A) -1 are compact as well. 

Since Ag is a bounded perturbation of Qg, it is immediate that Ag is closed. In 
order to show that Ag has compact resolvent, it suffices to show that (Ag — A) -1 
exists and is compact for some A G C (see, e.g., |Kat761 Theorems IV. 1.1 and 
III.6.29]). For A > 0, we can write 

(4.7) Ae-X=(l + n(Qg-Xy 1 )(Qg-X). 

By Lemma 13.51 and by Lemma 13.131 (i). (1), we have the uniform estimates 

A + Ai(0) 

llD^ + A)- 1 !! =\\DA 9 - 1 \\\\Ag(Ag + X)- 1 \\ < -==, 

v x n d ) 

Together with (|4.5p . it follows that Ag^ a (Ag + A) -1 is uniformly bounded for A > 0. 
Thus, p~6|) and the fact that \\(Ag + A) _1 | ->• for A ->• 00 and 9 G [0, 00] show 
that ||(Qe — A) _1 j| — > for A — > 00 and G [0,oo]. As a consequence, we can choose 
Ao > sufficiently large such that ||7?.(<2e — A) _1 || < 1 for A > Ao and so 

(Ae-xy 1 = {Qg+n-xy 1 = (Q e - xy 1 (i +n(Qg - xyy 1 

is compact for A > Ao. 

The symmetry of the (point) spectrum of Ag is evident since the entries of the 
block operator matrix Ag are differential operators with real coefficients. □ 

Although Ag is a bounded perturbation of Qg, the norm of the perturbation 
being ||a||, we cannot conclude that a(Ag) lies in a ||a|| -neighbourhood of a{Qg) — 
<t(— Ag) U a(— Aoo) since Qg is neither self-adjoint nor normal. Therefore, in the 
next proposition, we use a Neumann series argument to exclude points A G C from 
the (point) spectrum of Ag . 

Proposition 4.2. Let 9 G [0, 00] . If X ^ <r(-Ag) U a(- A^) and 

(4.8) \\a(A 00 + Xy\aA e ~ a'D)(A e + A)" 1 )) < 1, 
then X G p(Ag). 
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Proof. By IpTf) . we have A € p{Ae) if and only if A G p(Qe) = p(-^e) H /o(-Aoo) 
and J + lZ(Qg — A) -1 is boundcdly invertible. By (|4.6[) . the operator 

cA / -(Ae + A)- 1 

v Oj l v -(A 00 + A)- 1 A e , a (^ + A)- 1 -(A^ + A)- 1 

'7 - aiA^ + X)- 1 A et0l (A e + A)" 1 -a^ + A) -1 



I + n{Q 8 - A)" 1 =7 
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is boundedly invertible if and only if the operator J — o^-Aqo + A) 1 Ag^ a (Ag + A) 1 
is boundedly invertible. The latter holds if condition (|4.8p is satisfied. □ 



Remark 4.3. The block operator matrix can also be decomposed as 
Ag := 



-Ag 
-A c 



a 
Ag. a 



Here the first term is self-adjoint, while the second term is only relatively bounded 
with respect to the first one. A corresponding Neumann series argument leads to 
the same condition (|4.8[) for a point to be in the resolvent set. 



In the following we use Proposition ^. 21 to obtain an enclosure for the eigenvalues 
of Ag. To this end, we estimate the norm on the left hand side of (|4.8p by means of 
Lemma 13.71 According to the different resolvent estimates therein, we decompose 
the complex plane C as 

(4.9) C = Z U Z\ U Z 2 U Z 3 U Z 4 U Z 5 U Z & 

into the pairwise disjoint sets 



ZjQ .= 

Zi := 

z 2 ■■= 

Z 3 ■= 
Z 4 := 

Z 5 := 



z G C : Rez < -Ai(0), Imz = 



z G C : Rez < — Ai(oo), Imz ^ >, 



z g 



Z G 



Z G 



Z G 



oo, -Ai( 



- Ai(oo) < Rez < 0, Imz ^ 0, \z\ > A x (0) j, 

- A x (6») < Rez < 0, |^ + Ai(6>)/2| > A x (6>)/2, |r. ; \ y {H)\. 
-Xi(9) < Rez < 0, \z + Xi(9)/2\ < X x {6)/2 
Rez > 0, \z\ < Ai(6») 



z G C : Rez > 0, \z\ > Xi(6) 

The sets Z l are shown in FigureQ]for 6 = 1 = 1 where Ai (6) = tt 2 « 9.87, Ai (oo) w 20.19. 
Proposition 4.4. 7e£ # G [0, oo]. TTien 

(4.10) a(A) C E e := (-oo, -Ai(0)] U {A G C : /(A) > l} 
where the function f : C \ (— oo, — Ai(6 1 )] — > [0, oo) is defined as 

(4.11) /(A) := /j(A), X € Zj, j = 1,2, ... ,6, 
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Figure 1 . Decomposition (|4.9[) of the complex plane for I = 1. 



with 



/i(A) 
/a(A) 
/ 3 (A) 
/ 4 (A) 
/s(A) 
/e(A) 



|A| 



|ImA| 

|Aj 
|ImA| 

|A| 
|ImA| 



|A+|A||; |ImA|' 



"|A + |A||y |A + Ai(oo)| 



|A + A x (0)| j |A + Ai(oo)|' 

... ,„ VW) \ 1 



|A + Ai(0)| 



a\ \ a 



+ a a 



+ IHI IKI 



\X + Xi(e)\ J |A + Ai(oo)| 

|A + A x (0)| j |A + Ai(oo)|' 
1 



Ae4 
Ae4 

A e z 3 , 
, A e ^4, 
A e Z B , 
A e z 6 . 



|A + |A||y |A + Ai(oo)|' 
The set T,g is symmetric to M and bounded from the right with 
(4.12) a e := maxReS > -Ai(0), 

but it is not bounded from above nor from below in the complex plane C. 

Proof. By Lemma [d. 51 we have Ij-D^g ^ 2 || < 1- Hence condition f|4.8[) in Proposi- 
tion 14.21 is clearly satisfied for A ^ (— 00, — Ai(#)] if 



< 1. 



(4.13) IK^ + A)- 1 !! (HalPH^^ + A)- 1 !! + ||a|| ||a'|| WA^iAg + X)- 1 ^ 

Now we obtain the spectral inclusion (|4.10p from Proposition 14.21 by combining 
the various resolvent estimates in Lemma \3. 71 e.g., the estimates (ii), (b), and (2) 
for Z x . 
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The symmetry of Eg follows since Xi{8), Ai(oo) G R and hence /j(A) = /j(A), 
A ^ C\(— oo, —Xi(6)]. Because /e(A) — > for Re A — > oo, the set Eg is bounded from 
the right. The inequality (|4.12[) is immediate from the inclusion (—00, — Ai(0)l C Eg. 
If Eg were bounded from above or from below, there would exist an Mq > such 
that /i(A) < 1 for all A e C with |ImA| > Mo and Re A < — Ai(oo); on the other 
hand, 



flW - M JltxJ- M Mf^^ ReA^-cx;, 
a contradiction. □ 



In the sequel we describe the shape and the location of the set Eg in dependence 
of ||a|| and ||a'|| for a ^ 0; if a = 0, then clearly Eg = ( — 00, — Ai(6>)] . 

Lemma 4.5. Let f be given by (jlTTj) . Thentp: R 2 \((-oo, -Ai(0)]x{O}) -> [0,oo) 
defined by 

ipfori) :=/(A), e + i»? = AeC\(-oo,-A 1 (0)] ) 
is continuous and continuously differentiable on R 2 \ ((— 00, — Ai(#)] x {0}) with 

(4-14) !|(£, r/ )<0, ^(£,v) <0, CGR, ry>0, 



(/// 



d 



(4.15) ||(£,0)<0, |^(e,0) =0, f€R\(-c»,-Ai(fl)]. 

Proof. According to the definition of / in (|4.11[) . the function if is given by 

<P(£,V) = Wfori) •= A = e + i??e^, j = 1,2,..., 6. 

Hence ^ is continuous since the functions ipj are continuous on their domains Zj 
and, by the last claim in Lemma 13.71 the continuous extensions of two functions ifi 
and ifj coincide on common boundary points of Zi and Zj . 

It is elementary, but very tedious, to calculate the partial derivatives of the 
functions ipj on Zj and to show that the first partial derivatives of ifi and tpj coincide 
on common boundary points of Zi and Zj . This proves that (p has continuous first 
partial derivatives on R 2 \ ((—00, —\i(6)] x {0}) and is hence totally differentiable 
there. These calculations also show that (|4.14l) and (14. 15)) hold. □ 

In the next theorem the boundary of the set Eg is described by a function he] 
in particular, we derive formulas for the rightmost point ag of Eg (see (|4.12|l ). 

Theorem 4.6. Let 9 £ [0, 00] and a ^ 0. Then there exists a continuous strictly 
decreasing function hg: (—00, ag] —¥ [0, 00) with 

(4.16) lim hg(t) — 00, hg(ag) — 0, lim h' e (ag) = —00 

t— > — 00 t/*a s 

such that 

(4.17) o-(Ag) C Eg = {A G C : Re A < ag, |ImA| < «g(RcA)}. 
Depending on \\a\\ and \\ct'\\, the following cases occur: 
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(i) o<|H| 2 + HI|a/|1 



< Ai(oo): Then 



-Ai(6») < a g < 



and a$ is given by 

Ai(oo)+Ai(0) 



\ 



M°°)-M»h , Ai(9) | ||Qll2 



/••x n n9 ll a ll ll a 'l! » / \ , I, m IMI lla'l 
(ii) ||a|l + ILJULJi > Ai(oo) and ||a|| 2 + 



AM 
and is given by 



0<a e < Ai(6») 



AP)/ 

< Ai(oo) + Ai(6>): T/ien 



AjM + A^HHI 2 ^ AxM-Ax^)- ||a|| 2 V, / — - 

moreover, E# contains the disc Z4 — {z £ C : |z + Ai(#)/2| < Ai(#)/2}. 



(in) ||a|| 2 



> Ai(oo) + Ai(0): T/ien 



Ai(0) < a e 

and ag is the (unique) solution of the equation 
(4.18) 2%/A(A + Ai(oo) — ||a|| 2 ) = ||a|| ||o/| in (Xi(9),oo); 

moreover, E# contains the disc Z3 U Z4 U Z$ = {z € C : |z| < Ai(#)}. 

In Figure [2] the boundary of the set Eg containing the spectrum of the dynamo 
operator is displayed for the three cases above and 1 = 1, keeping the colour scheme 
for the sets Zi from Figure [TJ 



4(1 



20 



-20 



-40 















_J0_— ' ~~ 











-100 -80 -60 -40 -20 20 



Figure 2. Boundary of E e in Theorem 1461 (i). (ii), (iii) for I = 1 
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Remark 4.7. The second inequality in case (ii) can also be written as 

H 1. 2 , IHI IKII ^ , f T>,\ (a \ , IHl IKII 
a H < Ai(oo) + Ai(t/) H , 

This shows that case (ii) indeed appears and that (i), (ii), and (iii) exhaust all 
possible cases for ||a|| and ||o/||. 

Corollary 4.8. The operator Ag has no spectrum in the closed right half-plane if 

W 2 + P<A 1 (oo). 



Proof of Theorem 14.61 The existence of the function hg with the claimed 
properties follows from Proposition [431 Lemma 14.51 and the Implicit Function The- 
orem, applied to the restriction of if to the open upper half-plane. 

The monotonicity properties of ip in Lemma 14.51 induce corresponding mono- 
tonicity properties for the function /; in particular, / is strictly decreasing on 
(— Ai(0),oo). Therefore, since 

lim f(t) = lim f 4 {t) = oo, lim f(t) = lim / 6 (i) = 0, 

the equation /(A) = 1 has exactly one solution in (— Ai(0), oo); this solution is the 
unique zero ag of hg and hence equal to max Re Eg. The location of ag is classified 
by the three cases in Theorem 14.61 
(i) The condition in (i) is equivalent to /(0) = /4(C)) < 1 and thus ag £ (— Xi(0), 0]. 
Now the formula for ag is obtained by solving the quadratic equation 



(ag + X 1 (oo))(a e + X 1 (e)) = \\a\\ 2 Xi (9) + \\a\\ \\a'\\y/Xi(0) 
for ag in the interval (— \i(0),0], which is equivalent to fsiag) = !■ 

(ii) The first condition in (ii) is equivalent to /(0) = /5(C)) > 1, while the second 
condition in (ii) is equivalent to /s(Ai(0)) < 1. In this case ag £ (0, Xi(0)]. 
The formula for ag is obtained by solving the quadratic equation 



(ae + \x(oc)-\\af)(a e + \x(6)) = \\a\\ Wa'W^Me) 

for ag in the interval (0, Ai(#)], which is equivalent to fi(ag) = 1. 

Moreover, for A £ Z 4 = {z £ C : -Ai(0) < Re z < 0, \z + X 1 (9)/2)\ < \i(6)/2}, 

we have the estimates |A + Ai(0)| < Ai(0), |A + Ai(0)| < Ai(oo), and hence 




by the first condition in case (ii). This proves that Z4 C Eg. 
(iii) The condition in (iii) is equivalent to /(Ai(0)) = fe(\i(9)) > 1. Therefore 
ag £ (Ai(0),oo) and the equation /(A) = 1 is equivalent to /6(A) = 1, which 

is Km . 

Furthermore, for A £ Z 5 = {z £ C : Re z > 0,\z\ < Xi(9)}, we have the 
estimates |A + \i(8)\ < 2A X (0), |A + A x (0)| < Ai(oo) + A x (0), and hence 



./>,(A) > ( \\ar + \\a\\\\a'\\ i_ ) — / > 1, AeZ 5 , 

2y/\i(9) J Ai(oo) + Ai(0) 

by the condition in (iii). This shows that C Eg. From the monotonicity 
properties of / it follows that also Z 3 U Z4 c Eg. □ 
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Remark 4.9. The spectral inclusion in Theorem 14.61 does not provide a uniform 
bound for the imaginary parts of the eigenvalues in the left half-plane. Such a 
bound is obtained in the next section where we show that the imaginary parts of 
all eigenvalues are bounded by \\a'\\ (see Theorem 15. 10p . 

While Proposition 14.21 was applied in the proof of Theorem 14.61 to establish a 
global inclusion for all eigenvalues of Ag, it can also be used to study the local 
behaviour of eigenvalues. 

The eigenvalues of Ae for a = coincide with the eigenvalues of the diagonal 
elements — Ag and — Aoc and are hence real; their multiplicity is 1 for 8 G [0, oo) 
and 2 for 9 = oo. The following proposition provides a condition guaranteeing 
that, e.g., if 6 G [0, oo), for a ^ an eigenvalue remains on the real axis. With 
regard to the dynamo problem (jl.2[> . such a results could be called a local 

non-oscillation theorem. 

Proposition 4.10. Let 6 G [0, oo] and Xq G <t(— Ag) U a{— Aoo). Set 

S := i dist (X , (cr{-A e ) U <t{-A, x )) \ {A } 
and denote by Tq the circle centered at Xo with radius So. If 

(4-19) N| 2 + <TV ^r, 

y/\\o\+26 |Ao|+2V 

then, for 9 G [0,oo), the operator Ag has exactly one eigenvalue within the circle Tq, 
and this eigenvalue is simple and real; the operator A^o has exactly two eigenvalues 
within the circle Tq (counted with multiplicities), and these eigenvalues are either 
real or form a complex conjugate pair. 

Proof. For A G r , we have (see ([5715)) ) 



max \\(Ag + X) *\\ = — 
Aer d 



|A | +5 a 



max \\A e {A e + Xy L \\ = max X(Ag + X)- L \\ < 1 , 
Aero Aer oq 

max \\Al /2 (Ag + A) _1 | = max \\Ag(Ag + X)- 2 \\ 1/2 

^maxdK^+Ar^l + lAHK^+A)- 2 !!) 172 
Aei o 



< V|A |+2(5o 
So 

Using these estimates, it is not difficult to see that condition (|4.19j) guarantees that 
condition (|4.13l) . and hence condition (|4.8p in Proposition ^. 21 holds for all A G r . 
Now Proposition 14.21 yields that r C p(Ag). 

If we introduce the operators Ag := Qg + elZ for < e < 1, then the same 

arguments as above yield the inclusion r C p{A e e ^) for all < e < 1. Since the 
eigenprojection 

P(f\ = - 

2vri 



r 
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onto the part of the spectrum of Ag inside Tq depends continuously on e, the 
dimension of its range is constant for < e < 1 (see [Kat76[ Lemma 1.4.10, The- 
orem IV. 3. 16 and its proof]). For 9 e [0, oo), the operator A g , and hence every 
operator A g e ^ for < e < 1, has exactly one eigenvalue of multiplicity 1 inside Tq. 
The operator A^ , has exactly one eigenvalue of multiplicity 2 inside Tq and hence 
every operator A^ for < e < 1 has exactly two eigenvalues in To counted with 
multiplicities. 

Since all entries of the operator matrices Ag are differential operators with real 
coefficients, the spectra cr(Ag) are symmetric to R for all < e < 1 (compare The- 
orem l4.1[ ). Thus, for 9 € [0, oo), the single eigenvalue of Ag = Ae inside r must 

be real, while As) = Aoo has two eigenvalues inside Tq counted with multiplicities 
which are either real or form a complex conjugate pair. □ 

Corollary 4.11. Let 9 £ [0, oo) and let Fo be the circle around the largest eigenvalue 
— Xi{9) of the diagonal elements of Ae with radius <5o = (Ai(oo) — Ai(#)) /2. If 

2 HalllKH ^ (Ai(oo) - \i(9)) 2 



(4-20) n nnn < 

VAi(oo) 4Ai(oo) 
then the operator Ae has exactly one real eigenvalue in Tq. 



5. A SIMILARITY TRANSFORMATION AND A SECOND EIGENVALUE ESTIMATE 

Our second approach to study the spectral properties of the operator Ae is based 
on a quasi-similarity transformation of Ae . The transformed operator Be turns out 
to be a bounded perturbation of a self-adjoint operator Se- This allows us to prove 
another estimate for the eigenvalues of Ae which shows, in particular, that the 
imaginary parts of all eigenvalues have modulus at most ||a'||. In addition, we 
investigate the number of positive eigenvalues of the unperturbed operator Se- 

Proposition 5.1. Let9e [0, oo] and let the linear operator We in £2(0, 1)0^2(0, 1) 
be given by 

Wg := , V{W e ) := V{A X J 2 ) © L 2 (0, 1). 

Then the linear operator Be in £2(0, 1) © £2(0, 1) defined as 

(5.1) Bg := WgAgWg 1 

is closed and the eigenvalues of Be coincide with those of Ae, 

a p (B g ) = <j p (Ae); 

moreover, (rEj)*L is a Jordan chain of Ae at X E a v (Ae) if and only if (We Xj)j = o 
is a Jordan chain of Be at X. 

Proof. Since We is boundedly invertible and Ae is closed, the product Be is a closed 
operator (see |Kat76( Section III.5.2]). 

In order to prove the equality of the point spectra and the claim about the 
Jordan chains, we observe that, by definition, the domain of Be is given by 

(5.2) P(B e ) = {yeZ 2 (0,l)©Z 2 (0,l): W g 1 yEV(Ae), AgW^y eV(W e )} . 
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Let A € cr p (Ae) and let (xj)j =0 C T> (Ag) be a corresponding Jordan chain, that is, 
(Ag — X)xj = Xj-i for j = 0, 1, . . . , k, where x_i := 0. Since V(Ag) C V(Wg), it 
follows that Agxj = XqXj + Xj_i <E T>(Wg). Hence Wex.,- € V(Bg) by ([53]) and 

(B fl - \)WgXj = (WgAgWg 1 - \)WgXj = (WgAg - XWg)x 3 = WgXj-x 

for j — 0,1,..., k. Vice versa, let A £ <J P {Bg) and let (yj)| =0 C T>(Bg) be a 
corresponding Jordan chain, that is, {Bg — X)yj = yj-i for j = 0, 1, . . . , k, where 
y-l := 0. Since W^Vj &T>(Ag) by ([Oft , we obtain 

(A - Ajw- 1 ^ = (>%-%>v e - Ajyiv 1 ^ = (nv% - AW e ->, = >%- Vi 

for j = 0, 1, . . - , Jfe. □ 

In the next theorem we show that the operator Bg is a bounded perturbation of a 
self-adjoint semi-bounded block operator matrix. First we need an auxiliary lemma. 

Lemma 5.2. Let 6 e [0, oo]. Then 

1 /2 

i) a A g is Ag- compact, 

1 /2 

ii) Ag a is Aoo-compact. 

Proof, i) Since Ag has compact resolvent by Proposition 13.11 ii) , we conclude that 
aA^Ag 1 = aAg 1 ^ 2 is compact, 
ii) First we show that 

(5.3) V(A l J 2 )<zV(A 1 J 2 a). 

The description of V(Ag^ 2 ) in Proposition 13. H i) and the assumption a€C 1 ([0,l]) 

show that xeV(A 1 g /2 ) implies ax€V(A 1 J 2 ), that is, V(A 1 J 2 ) C V(A 1 g /2 a). Since 

V{A]i 2 ) = {x e V(A X J 2 ) : x(l) = 0} C V{A X J 2 ) by PropositionEUi), the inclusion 
dOt follows. 

-j^ / 2 1/2 1/2 1/2 

Because Aq£ and .Ag a are closed, (|5.3p implies that Ag aAoo is a bounded 
operator (see |Kat76[ Remark IV. 1.5]). Since A^ has compact resolvent by Propo- 
sition [3TT] ii) , it follows that A^ 2 aA^ = J^J 2 aA^^ 2 A^^ 2 is compact. □ 



Theorem 5.3. Let 6 € [0, oo]. Then the operator Bg defined in (|5.ip has the form 

V{B e ) =V(A ) ®V(A oa ). 



Ag Al /2 a 
A' 1 ' 2 




V(Sg) :=V{A e )@V{A 00 ), 

V{Tg) :=L 2 (0,1)®L 2 (0,1), 

thenSg is self-adjoint and bounded from above with compact resolvent, Tg is bounded 
with ||7g|| < \\ a '\\j an d 

(5.7) Bg = Sg + Tg. 
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Proof. First we prove the identity (|5.4|l . By the definition of Ag, a in Proposi- 
tion EH we have V{A e , a ) = V(Ag) and thus V(A e ) C £>(Ag /2 ) = V{A e . a A~ 1/2 ); 

by Lemma [5.21 ii) (see also (|5.3p ). we have T>(Aoq) C V^A^J 2 ct) . This shows that 
T>(Ag) © 25(^00) is the domain of the block operator matrix in (|5.4I) . Formally, the 
relation (|5.4j) follows from the definition of in (|5.1|) ; it remains to be shown that 
V{B e )=V(Ag)(BV(A 00 ). 

According to (15. 2|) . we have y £ X>(Be) if and only if W e 1 y £ V(Ag) and 
AeW^y £ V(W e ). Since D(^) = V(A e )®V(A 00 ) by ([221) and 

(5.8) A e , a A e 1/2 = uA x J 2 - a'DA 8 1/2 

by ([3~T2"|) . it follows that y = (y 1 ,y 2 ) t £ V(B g ) if and only if 
f (a) G D(A e ), y 2 £ X>(4»), 

\ (b) -^/ 2 yi + ay 2 el?(^ /2 ), (a^ /2 - a'£>A- 1/2 ) yi - A M y 2 £ L 2 (0, 1). 

1 /2 

The first condition in (a) is equivalent to y± G *D{AJ ). Thus the second condition 
in (b) always holds as DA g is bounded by Lemma T3. 51 By Lemma T5.2I ii) (see 

1 /2 

also (|5.3p ). we have T>{A ca ) C T>(A g a) and hence the first condition in (b) reduces 

to A\ /2 yi £ V(Al /2 ). Altogether, y = {y 1: y 2 ) t £ V{Bg) if and only if Vl £ V{A g ) 
and y 2 £ V{A oa ). This completes the proof of (|5.4[) . 

To see that Sg is symmetric, we observe that a £ C([0, 1]) and a is real-valued. 
This implies that the corresponding multiplication operator is bounded and self- 
adjoint and hence 



(a<y = K /2 )v=4 l/2 



Furthermore, Lemma 15.21 implies that 
So = 



a. 



A g \ . ( Ai /2 a 



-^r^ 3 y 

is a relatively compact perturbation of its block diagonal part diag(— Ag, — Aoo). 
The latter is self-adjoint and semi-bounded (in fact, negative) with compact resol- 
vent by Proposition 13.11 Hence Sg has the same properties by the Rellich-Kato 
Theorem, the stability theorem for semi-boundedness, and Weyl's essential spec- 
trum theorem (see |Kat761 Theorems V.4.3 and V.4.11], |RS781 Corollary XHI.4.2]). 

Finally, because a' £ C([0, 1]) and DAg 1 ^ 2 is bounded with HZM^ 1 ^ < 1 (see 
Lemma I3~5|) . the operator Tg is bounded with ||7e|| < 

Since V(Sg) = V(Bg), the identity dSTTJ) is immediate from (JOJ) and (fSTSj) . □ 

Remark 5.4. The boundedness from above of Sg is proved independently in Propo- 
sition [^H where also a concrete upper bound for Sg is established. 

Proposition 5.5. Let 9 £ [0, 00]. The spectrum o~(Sg) = a p (Sg) of the self-adjoint 
operator Sg satisfies 

(5.9) a p {S e )\a p (-Ag)=a p (S 2 ,g) 

where the so-called Schur complement S 2 _g is the operator function given by 

S2, (X) := -A x - A + aAg /2 (A e + Xy 1 A 1 J 2 a, V(S 2 , e {\)) := V(A 00 ), 

for all X £ C\ o-{—Ag). Moreover, the multiplicities of the eigenvalues of Sg and 
of S%fi coincide. 
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Proof. Let A ^ a(—Ag) = cr p (— Ae). It is not difficult to check that the block 
operator matrix Sg can be factorised as (see, e.g., |Tre081 Theorem 2.2.14]) 

(5 - 10) S9 - X -{-aAy 2 (A e +\)-i l)[ S 2 ,eWA0 I 

Since the outer two factors are bounded and boundedly invertible and so is — Aq— A 
by the assumption on A, it follows that Sg — A is boundedly invertible or injective 
if and only if so is the operator S 2 ,g(\) (see, e.g., |Tre08[ Theorem 2.3.3]). This 
proves (|5.9[) . 

It is easy to see that y = (2/0,1,2/0,2)* & D(Sg) is an eigenvector of Sg at 
•^0 £ a p (Sg) if and only if 2/0,2 € 2?(Aoo) is an eigenvector of S 2t g at Ao, that 

is, 52,0^0)2/0.2 = 0, and 2/0,1 = (A# + A)~ 1 Ay 2 a 2/0,2- Clearly, the self-adjoint op- 
erator Sg has no associated vectors. Since S2 is a self-adjoint operator function, 
that is, S2,g(X) = S 2l g(X)* for A e C \ a(—Ag), and its derivative 

S' 2 ,e( X ) = - J - a^ /2 (A e + X)- 2 Al /2 a < -I 

is uniformly negative for all A G C \ cr(— Ag), S^g has no associated vectors either 
(compare !Mar881 Lemma 30.13]). □ 

Next we establish an explicit upper bound for the self- adjoint operator Sg. To 
this end, we first show that Sg satisfies the following two-sided estimate, which 
might be of independent interest. 

Lemma 5.6. Let 9 £ [0,oo] and let 7 S (0, 1] be arbitrary. Then the operator Sg 
admits the two-sided estimate 



\ /-I-7U0 




(5.11) 1 „ < Sg < 1 o 

-A x a 2 \ -A^-a 2 



7 / \ 7 



in particular, for 7=1, 



f-2A e \ /0 

(512) 1 -A.-*)**^ 

Proof. Let y = (2/1,2/2)* & V{Sg) = V(Ag) 8 £>(^oo)- Then the right inequality in 
(15.111) follows from the estimate 

V2, 



(Sgy,y) = -(A e yi,yi) - (A 00 y 2 ,y2) + 2Re(A g ' a 2/2,2/1) 



< -(A e 2/i,2/i) - (^2/2,2/2) +27 1 / 2 ||^ /2 2/ 1 ||-i 



1/2 



I "2/2 I 



< -(Aeyi,yi) - {A QO y 2 ,y 2 ) +j{A e y 1 ,y 1 ) + -(a 2 y 2 ,y 2 ); 

7 

the left inequality in (15.11[) is obtained analogously. □ 
Proposition 5.7. If a is constant, a = ag, we also have the estimate 

(5.13) maxer^e) = maxa(Sg) < — 7— • 
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Proof. Lcty = (y u y 2 f € V(Sg) = V(Ag)®V(A 00 ); then y 2 € £>(AxO CP(^ /2 ) C 
V{A 1 J 2 ) by Proposition 13 . 1 1 i) . Since a is constant, it commutes with A^ 2 . Hence 

(S e y, y) = - {A e yi,yi) - (A QO y 2l y2) + (aoAl /2 y 2l yi) + (a Q Al /2 y l7 y 2 ) 

y2 ) + |ao|||^ /2 y2 |||| J/1 || 



l a oln ,,\ . |ao| 



< - ( ^(A e yi,yx) -^% 2 || ) + ^||z/ 2 || 2 



2 



/ — v m ii^y a ifai ii „\ , i«oi 2 iK^rf n ii ; 



<^(IMI s + llrf) = ^ll!/ll 2 

if we show that 

Il4 /2 y2 ||<||^ 2 y2 ||, y 2 eV(A oa ), 

or, equivalently, || Ag^ 2 A^^ 2 \\ < 1. This follows from the relation A X J 2 A^l 2 = 
(A^ 1/2 Al /2 )* (note that A X J 2 A~^ j2 is bounded as V{A\L 2 ) C X>(^ /2 )) and from 
the estimate 

IjA- 1 /^ 2 *!! 2 = {A^Ay 2 x,Al /2 x) < {A^A l J\,A x /\) = \\x\\ 2 , x € 2?(^ /2 ) ; 

where we have used the inequality A^ < Ag 1 from Proposition 13 . 1 1 iv) . 

Since a is constant and hence a' = 0, we have Bg = Sg and so (|5.13p follows 
from Proposition 15. II □ 

The following proposition provides an upper bound for the unperturbed operator 
Sg for arbitrary functions a. 

Proposition 5.8. Let G [0, co]. Then the self-adjoint operator Sg defined in (|5.6|) 
satisfies the estimate 

max<r( l Se) < sg 

where sg G [— Ai(0), — Ai(oo) + ||a|| 2 ] is given by 



so 



Ai(oo) + \i(9) , / / Ai(oo) - A 1 (6') N 2 



[-ArM + Hall 2 ~ z/ H 2 >Ax(oo). 



Remark 5.9. i) The bound sg satisfies 

-X 1 {9)<sg<0 \\a\\ 2 < Ai(oo), 

< s e < -Ai(oo) + ||a|| 2 <*=^> ||a|| 2 > Ai(oo). 

ii) If ||a|| 2 < Ai(oo), then sg can be written equivalently as 



sg 



-Aitf) + Hx/AiW tan flarctau^l^p^ 

V2 Ai(oo)-Ai(0) 
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iii) In the particular case 9 = oo, the expression for Soo simplifies to 



-AiM + IHIVAiM if ||a|| 2 <Ai(oo), 
-Ai(oo) + ||a|| 2 if ||a|| 2 > Ai(oo). 



= — Ai(oo) + ||o!|| max{||a||, y/ Ai(oo) }. 

Proof of Proposition 15.81 Since — Ag < — Ai(0) and — Aoo < — Ai(oo), the right 
inequality in (15.111) yields that, for arbitrary 7 £ (0, 1], 

(S e y,y) < -(1 - 7) Ai(0) ||yi|| 2 + ( - Ai(oo) + ±N| 2 ) \\y 2 \\ 2 

< max {^(7), ^(7)} l|y|| 2 

where we have set 

hi (7) :=-(l- 7 )Ai(0), h 2 (7) :=-Ai(c») + -||a|| 2 , 7 €(0,1]. 

7 

If ||a|| 2 > Ai(oo), it is not difficult to see that ^2(7) > > ^1(7) for all 7 £ (0, 1]; 
in this case, the optimal estimate is obtained for 7 = 1, that is, 

(Sey,y) < (-A 1 ( o) + ||a|| 2 )||y|| 2 . 

If ||a|| 2 < Ai(oo), a short calculation shows that the function maxj/i!^), /i 2 (7)} 
attains its minimum at the point 70 £ (0, 1] where hi and h 2 intersect and which 
is given by the relation 



1 / Ai(oo)-Ai(g) f X 1 (^) - X 1 (9) \ 2 

l0 = xA6){ 2 + Vl 2 J +M*)IN 

in this case, the optimal estimate becomes 



Ai(oo)+Ai(fl) , //A 1 (oo)-A 1 (0) N 2 



5 e y,y < - ±v ' ' 1V y + W 1V ' iy ' + Xi(0)H\ 2 \\y 



2 



where we have used that Ai (6) < Ai (00) for all 8 £ [0, 00] by Proposition ^. ll iii). □ 



While the first spectral enclosure in Theorem 14.61 only provides some region in 
the complex plane containing the eigenvalues, the next theorem gives more detailed 
information in terms of the eigenvalues of the self-adjoint operator Sg. 

Theorem 5.10. Let 9 £ [0, 00]. The eigenvalues of Ae lie in discs of radius \\a'\\ 
around the eigenvalues of the self-adjoint operator Sg: 

(5.14) a{A e ) C {A £ C : dist (A, a p (Sg)) < \\a'\\}. 
As a consequence, 

(5.15) a{A e ) C {AeC : ReA < ag, |ImA| < \\a'\\} U {AeC : |A - s \ < \\a'\\} 
and 

(5.16) \£a(Ae) => |ImA| < ||a'||, Re A < s e + \\a'\\ =: bg. 

Proof. All claims follow from Theorem [53] and Proposition l5.8l bv means of classical 
perturbation theorems for self-adjoint operators (see, e.g., |Kat76, Theorem V.4.5]) 
if we observe that Bg = Sg + Tg and ||7e|| < \\c/\\. □ 
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The following corollary guarantees that Ae has no eigenvalues in the closed right 
half-plane; it is an immediate consequence of Theorem 15 . 101 and of the definition of 
se in Proposition 15.81 

Corollary 5.11. Let 9 6 [0, oo]. Then A$ has no spectrum in the closed right 
half-plane if 

., „ n— — 7 „ „, Ai(oo) + Ai(0) / ( Ai(oo) — \ 2 , ,„.„ ..„ 
||a|| < VAiM, \\ol\\ < ' 2 ^ - J f > 2 + Xi(0)\H\ 2 - 

In a similar way as Proposition 14.101 the following local result can be proved. 

Proposition 5.12. Let 9 £ [0, oo] and let Ao G o~(S$) be an eigenvalue of Sg with 
multiplicity mo- Set So :— dist (Ao, cr(Sg) \ {Ao})/2 and denote by To the circle 
around Ao with radius So- If \\a'\\ < So, then the operator Ag has mo eigenvalues 
inside To (counted with multiplicities); if mo = 1, then the unique eigenvalue of Ae 
in To is real. 

While in Section [5] the spectrum of the unperturbed operator Qe is known, this 
is not true for the unperturbed operator Se. Therefore, in the remainder of this 
section, we investigate the spectrum of the semi-bounded self-adjoint operator Se 
by means of variational principles. 

The quadratic form version of the min-max characterization of the eigenvalues 
of self-adjoint operators bounded from below (see |WS721 Theorem 1, Section 2.2], 
[RS78| Theorem XIII.2]) shows that the eigenvalues X k (Sg), k € N, of S e enumer- 
ated such that X\(Se) > X 2 (Sg) > • ■ • and counted with multiplicities, are given by 

(5.17) Xj(Sg) = max mm s [y], j = l,2,.... 

CCQ(Sq) y€C 
dim£=j II y ||=1 

Here Q(Sg) = V(A l J' 2 ) V(AH 2 ) is the form domain of Sg and Se is the corre- 
sponding quadratic form given by 

(5.18) Sg [y] := - (A 1 /' l Vl , A 1 /" 'm) + (ay 2 , A l J 2 yi ) + (A l / 2 yi ,ay 2 ) - {A^y^A^) 

for y = (yi,y 2 ) t e Q(Sg) (see (JSHJ) and [RS80l Section VIII.6]). 

Proposition 5.13. The number ko of positive eigenvalues of the operator Se in 
1/2(0, 1) © £2(0, 1) coincides with the number of positive eigenvalues of the shifted 
Bessel differential operator — A^ + a 2 in £2(0, 1), and we have 

(5.19) 0<X J (Sg)<X J {-A oo +a 2 ), j = 1, 2, . . . , k . 

Proof. Denote by Kq the number of positive eigenvalues of —A^ +a 2 . The inequal- 
ity ko < Ko a n d the inequalities (|5.19[) are immediate consequences of the right 
estimate in (|5.12l) . which implies that 

Sg[y] = (Sgy,y) < ((-A 00 +a 2 )y 2 ,y2) 

for y = (2/1,1/2)' € T){Sg) = T>(Ag) V(A 00 ) and of the variational principle (|5.17p . 

It remains to be shown that fco > to- By the definition of ko and by the varia- 
tional principle (|5.17|) applied to the semi-bounded operator —Aoo+a 2 , there exists 
1/2 

a subspace C C V{AJ ) = Q(A oa ) (the form domain of the positive operator Aoo) 
such that dim£ = no and 

min ( - (A)l 2 x, A^x) + (ax, ax)) > 0. 

ll*ll=i 
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Using (|5.18[) . it is easy to see that the subspace j (A g 1 ^ 2 a 2/2,2/2)* : 2/2 € c\ C Q(<S#) 
has dimension k and 



(5.20) 5 6 



A e ay 2 



= ((-Aoc + a 2 )2/ 2 ,2/ 2 ) > 0, y 2 et(ZV(A 



1/2N 



This and the variational principle (|5.17[) show that X Ko {Sg) > and hence Sg has 
at least kq positive eigenvalues. □ 

The eigenvalues of the block operator matrix Sg coincide with the eigenvalues of 
its Schur complements 

Si,g(X) := -A e -X + Aj /a aC4» + X^aA^ 2 , V(S h g(X)) = V(Ag), 

S 2< g{X) ■= ~ A °°- A + aAl /2 {A e + Xy'A^a, V(S 2 ,g{X)) = Z?(4»), 

which are defined for A £ C \ cr(— Aqo) and A € C\ a(—Ag), respectively: 

a p (Sg) \ a^-Aoo) = er p (Si, e ), a p {Sg) \ a {-As) = cr p (S 2: g) 

(compare Proposition 15 . 5p . This allows us to characterize and estimate the eigen- 
values of Sg in the intervals (— Ai(oo), 00) and (— Ai(0), 00) by variational principles 
for S\_g and 62,0, respectively (see |EL04] and, e.g., |Tre081 Section 2.10]). As an 
example, we consider the eigenvalues of Sg in (— Ai(0),oo). 

Lemma 5.14. Let 9 £ [0, 00]. The Schur complement S 2t g satisfies the estimates 

62,0 (A) < -4, + A A c? - X, Ae (-Ai(0),O], 

S 2 ,g{X) K-A^ + a 2 - X, A e [0, 00). 

Moreover, the derivative of S 2 ,g is strictly negative and satisfies 

S' 2 ,g(X)<-I, Xe{-X 1 {6),oo). 
Proof. Let A € ( — Ai(0),oo). The Schur complement S 2: g can be rewritten as 

S2,e(A) = -Aoo - A + aAg(A e + X^a 

= -Aoo - X + a 2 - Xa(A g + A) _1 a. 

Since Ag > Xi{6) > (see Proposition 13. lj) . the resolvent of Ag satisfies the two- 
sided inequality < {Ag + A)" 1 < l/(Ai(0) + A). For A e ( — Ai(0),O] , the right 
inequality yields the first estimate claimed for S 2 ,g(X), while for A € [0, 00) the left 
inequality yields the second one. 

The inequality for the derivative of S 2t ${X) follows from the identity 

S' 2 g{X) = -I-aA 1 J 2 {Ag + X)- 2 A 1 e /2 a, X e ( - Ai(0), 00). □ 



Proposition 5.15. Let 9 € [0,oo]. For every subinterval [a, 00) C ( — Ai(#),oo) 
let Xk a (Sg) < ••• < Ai(5g) 6e eigenvalues of Sg in [a, 00) (counted with 
multiplicities) . Then 

(5.21) Xj{Se) = min max £2,0(2/2)) j = 1, 2, . . . , fc Q , 

dim £ = 3 || y 2 11=1 
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where P2, 0(1)2) is the (unique) zero of (52,e(-)z/2, 2/2) on [a, 00) if a zero exists and 
£•2,0(2/2) := —00 otherwise for 1)2 € V(A 00 ). Moreover, the number k a of eigenvalues 
of Sg in [a, 00) is given by 

k a = dim£[ 0iOO )(S , 2,e(a)); 

here Ci(S2,s(X)) denotes the spectral subspace of the self-adjoint operator 52,e(A) 
corresponding to an interval I C ( — Ai(0),oo). 

Proof. The Schur complement 52,0 satisfies the assumptions of |EL04[ Theorem 2.1] 
on every subinterval of ( — Ai(0),oo): 2?(52,g(A)) = V(A ao ) is independent of A, 
S2.6 is strictly decreasing with S' 2 g < —I, and (5 2j 6>(A)2/2, 2/2) — > —00, A — » 00, 
for all 7/2 € I^Aqo) by Lemma T5. 141 Now all claims follow from the fact that the 
eigenvalues of Sg in ( — Ai(0),oo) coincide with those of £2,8 by Proposition 15.51 
and from the variational principle in [EL041 Theorem 2.1] applied to 52, e- □ 

Remark 5.16. Setting a — in Proposition 15 . 1 51 we obtain another proof of Propo- 
sition [5331 In fact, the estimate (I5.19[) follows from the variational principle (|5.21[) 
if we observe that 52,6»(A) < — A x + a 2 — A for A e [0, 00) by Lemma f5. 141 and hence 

, , . {(-A 00 +a 2 )y 2 ,y 2 ) 
P2.e(V2) < -, \ -, V2 € V(A 00 ), 

(2/2,2/2) 

the right hand side being the zero of the function A 1— > ((— A, x + a 2 — A) 2/2, 2/2)- 
Further, since 5 2j e(0) = — Aao + a 2 , Proposition 15. 151 shows that the number k of 
positive eigenvalues of Sg is given by 

k = dim£[ 0:OO )(5 2 ,e(0)) = dim£[ 0iOO )(-^oo + a 2 ). 
6. Comparison of the two eigenvalue estimates 

The two eigenvalue estimates obtained in Theorem l4.6l and Theorem l5.10l show. in 
particular, that every eigenvalue of the dynamo operator Ag satisfies the inequalities 

|ImA| < ||a'||, Re A < mm{a e ,bg}; 

here the uniform bound for the imaginary parts was proved in Theorem 15 . 1 01 and ag 
and bg = sg + \\a'\\ are the right bounds for a(Ag) derived in Theorems 14. 6l and l5.10[ 
respectively. In this section we prove that, apart from a small bounded set of values 
of ||a|| and \\a'\\, we always have ag < bg. Hence, in general, a combination of the 
two eigenvalue estimates from Theorems 14.61 and 15.101 yields the best result. 

First we present two auxiliary technical lemmas. The graphs of the functions 
introduced therein are displayed in Figure [3] below. 

Lemma 6.1. The functions 



1 rn m-v m n 1 n\ If Ai(°°) - Ai(#) \ 2 
fci:[0,VAi(oo)]->[0,oo), h(t) := '- — -Wf '— — 1 + Xi(9)t 2 , 

fc 2 :(0, v / ATM]^[0,oo), k 2 (t) := X ^°°}~ ^ y/X^ff), 



Ai(oo)-i 2 + Xi(9) 



fc 3 :(0,v^M + A^)]^[0,oo), k 3 (t):= n ; L -^2^(9), 

are continuous and strictly decreasing with 

k(VMoo)) = MVM°°)) = fc 3 (AR+M9)) = 0. 
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They satisfy the inequalities 

fci < fc 2 on (0, y / X 1 (oo)), 

(6.1) 



fc 2 < fc 3 on (O lV /Ai(oo) + Ai(6»)), 



where we have set fe(t) := /or i € (-y/ Ai(oo), ^/Ai(oo) + Ai(0)]. Moreover, 

IKII < h(\\a\\), 



bg<0 
ag<0 

<a g < Ai(6») 



H < AM. 

||a|| < VAi(oo), 



l«'ll <*2(l|a||), 



||a|| < a/Ai(oo) + Ai(0), Ho'II < fcs(||a||). 



Proof. The claims for the functions fci, &2, and fc 3 are easy to check. The last 
two equivalences follow if we use the formulas for ag in Theorem 14.61 and solve 
the corresponding inequalities on the left hand sides for \\a'\\. The condition bg — 
sg + 1| a' || < is satisfied if and only if sg < and ||a'|| < —sg. It remains to use the 
respective formula for sg from Proposition 15.81 to obtain the first equivalence. □ 



Lemma 6.2. Define two functions k 4 implicitly by the equation 



for t £ [0, -^/Ai(oo)] and /ei 
fe 5 



V Ai(oo), + \/Ai(oo) - 



0,oo 



2 V v ' 4 

h{t) ■= A 1 (oo)-A 1 (0)+t v ^(0)-t 2 . 
T7ien t/ie graphs of k^, kj, and h form a continuous curve T ex connecting the points 
d := (VaTM.o), C 2 := (V^M, v/aTMATM-AiW), 



c 3 



(^ + /a iM -^o 



Le£ A ex 6e i/ie open bounded set surrounded by the curve F ex and £/ie segment C1C3 
on the ordinate axis. Then 



A ox c {(HJa'H) : fc 2 (||a||) < ||a'|| < fc 3 (||a||), ||a'|| < VAi(oo)Ai(0) - \i(6)} 

and, for < ag < Xi{9), we have 

bg = ag ^ (||a||,||a'||) G T cx , 

b e <ag ^ (H,||a'||) € A cx . 

Proof. It is easy to sec that the graph of the function £5, which lies to the right of 
the vertical line ||a|| = y/ Ai(oo), is strictly decreasing from the point C 2 on this 
line to C3, the zero of h. Moreover, if < ag < Xx(6) and ||a|| > \J Ai(oo), the 
formulas for ag in Theorem 14.61 (ii) and bg — sg + \\a'\\ = — Ai(oo) + ||a|| 2 + ||a'|| 
(see Proposition 15. 8[) show that h > on the interval where it is defined, 

bg < ag <^> \\a'\\ < k 5 (\\a\\), 

and that equality holds on the graph of h. 
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By elementary calculations (see Remark[6?3]below) , one can show that the graphs 
of the implicitly defined functions fcj and k^ form an arc from C\ to C2 lying to 
the left of the vertical line ||a|| = -\/Ai(oo). Furthermore, if < ag < Xi(9) and 
\\<x\\ < \/Ai(oo), the formulas for ag in Theorem l4.6l (ii) and for bg = sg + \\a'\\ with 
sg as in Proposition 15.81 show that 

bg<ag ^ kl (H) < \\a'\\ < kt(\\a\\) 

and that equality holds on the graphs of kf . □ 

Remark 6.3. The two functions fc 4 in Lemma 16.21 can be calculated explicitly: 



^(t)=^+^v^w-i/( A ' M : A ' ( " ) ) 2 +A ll «i, i 



for t e [0, ^/Ai(oo)] such that the term under the last square root is positive. The 
set of these t is a (small) interval of the form [/x, -\/Ai(oo)]. A formula for fi may 
be found, e.g., with MAPLE, but it is extremely involved; a lower bound for /1 is 
the point where ki attains the value sj Ai(#)Ai(oo) — Ai(#) (the height of A ex ): 



(6.2) „ > _^M±^m + ^VhM+vmy + AlM . 

Proposition 6.4. Lei G [0,oo]. XTien i/ie rig/ii bounds ag and bg for cr(Ag) 
established in Theorems 14.61 and \5.Wl respectively, satisfy 

a 6 <be ^ (|| a ||,|| a '||) e ([0,oo) x [0,oo)) \A ex , 

a e = bg (||a||, ||a'||) G T cx , 

where the bounded set A ox and the curve T ex are as in Lemma \6.2[ Moreover, 

a e <0<b e ||a|| 2 < Ai(oo), fci(||a||) < ||a'|| < fc 2 (||a||); 

the last condition can be written equivalently as 



M°°HM*) _^ Mo.)-A.w y + ^ w|a| , < ||Ql| < Ms^jjM!!^. 

Figure [3] illustrates the sets occurring in Proposition ^. 41 the exceptional set A ex 
where bg < ag is the highlighted small bounded set between the graphs of ki and 
k%\ the set where ag < < bg is the highlighted unbounded set enclosed by the 
graphs of ki and fc 2 . 



Remark 6.5. Proposition 16.41 shows, in particular, that we always have ag < bg if 
one of the following holds: 



i) IKH > ^Ai(oo)Ai(6») - Ai(6>), 



ii) \\a\\ < - 



...s II n ^ y/W) , I \ ( \ 3Ai(g) ' 
m ) INI > o + V Ai(oo) — . 
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This follows from Lemma l6.2l if we observe that the lower bound in i) is the maximal 
value of ||a'|| in A ex , the upper bound in ii) is the lower bound for the left end-point 
/j of the domain of definition of the functions kf (see Remark l6.3|) . and the lower 
bound in iii) is the right end-point of the domain of definition of the function k§ . 

Proof of Proposition 16.41 In order to compare ag and bg, we have to dis- 
tinguish the cases (i), (ii), and (iii) in Theorem 14.61 for ag and the two cases 
IMI < \/Ai(oo) and ||a|| > \J Ai(oo) for bg = sg + \\a'\\ according to the defi- 
nition of sg in Proposition 15.81 

Elementary but lengthy and tedious calculations show that in the cases ag < 
(Theorem 14.61 (i)) and ag > Xi(0) (Theorem 14.61 (iii)). the equation ag = bg has no 
solution and ag < bg. Moreover, Lemma f6 . 1 l implies that the cases ag < and bg > 
appear simultaneously if and only if || a| | < ^/Ai(oo) and fci(||o;||) < \\a'\\ < fc2(|M|) 
since we have fci < k 2 by (|6.ip . 

In the case < ag < \i(6) (Theorem 14.61 (ii)). Lemma T6.2I shows that ag = bg 
on the curve T cx and ag < bg if and only if (||a||, \\a'\\) A cx . □ 

In Figure [4] the boundaries of the two spectral enclosures from Theorem 14.61 
and Theorem 15.101 are displayed together; the set from Theorem 14.61 is the set 
with unbounded imaginary part and smaller real part < ag, the enclosing set from 
Theorem 15.101 is the set with bounded imaginary part and larger real part < bg. 
Here the values ||a|| = 1.5 and \\a'\\ = 15 are chosen such that ag < < bg; the 
corresponding point (||a||, ||a'||) = (1-5, 15) is marked in Figure [3] by a black dot. 

Proposition l6.4l shows that, in general, the upper bound of Theorem 14 . 61 improves 
the estimate of Theorem 15.10) only for values (||a||, ||a'||) £ A ex , the enclosure of 
Theorem [540] is better. The following corollary is a direct consequence of all these 
results. 

Corollary 6.6. Let 9 £ [0,oo], let ag, bg be the bounds established in Theorems 14.61 
and 15.101 respectively, and let A cx C [0, oo) x [0, oo) be the bounded set defined 
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2(1 



INN1.S, ||a'||=15 




-20 



20 



Figure 4 . Spectral enclosures of Theorems 14.61 and 15.101 (ag < < bg) . 



in Lemma 16.21 Then every eigenvalue A of the dynamo operator Ag satisfies 



|ImA| < 



Re A < 



if (IMIJo'll) e ([0,oo)x [0,cx))) \A C 
if (||a||,||a'||) 



The case that the bound for the real part is less than is of particular inter- 
est from the physical point of view since the onset of the dynamo effect requires 
supercritical modes, that is, eigenvalues with real part greater than 0. 

Proposition l6.4l shows that, if ag < 0, then we always have ag < bg. In particular, 
this shows that Corollary 14. 81 is stronger than Corollary 15.111 for the special case 
6 = I this gives: 

Anti-dynamo theorem. The dynamo operator Ai has no spectrum in the closed 
right half-plane if 



(6.3) 



IN 



< Ai(oo). 



Note that, since the above result has been obtained by operator theoretic esti- 
mates, a violation of the condition (16.31) provides only a necessary, not a sufficient, 
condition for the existence of supercritical dynamo regimes. 

Remark 6.7. By Lemma 13.31 we know that \i(l) and Ai(oo) are the first non-zero 
zeros of the Bessel functions J/_i/2 (VA) and </;+i/2 (V^A), respectively. If we indicate 
this dependence on the parameter I £ N for a moment and write Ai (Z) = Ai (I; I) and 
Ai(oo) = Ai(oo, I), we see that Ai(/i; h) < Xi(J,2', h) and Ai(oo; h) < Xi (oo; I2) for 
h,h € N with h < l 2 (see [AS641 9.5.2]). Therefore, if condition ((673)) is satisfied 
for l\ G N, then it is also satisfied for I2 G N with I2 > l\. 

This provides an explanation of the numerical observation that, for increasing 
|ja|| and \\o/\\, criticality usually starts from dipole modes (1 = 1) and quadrupole 
and higher-degree modes (I > 1) become supercritical only for larger values of ||a|| 
and Ha'H (see [SG03] 1. 



SPECTRUM OF MAGNETOHYDRODYNAMIC MEAN-FIELD DYNAMO OPERATORS 31 



7. Examples 

In this section we illustrate our results by applying them to the physical dynamo 
problem (jl.ip . (|1.2p and to the idealized dynamo problem (jl.lj) . (| 1 . 3|) . Here we 
have to set = I and — oo, respectively, in the previous statements. 

The simplest case of constant a, where all eigenvalues are real, already shows 
that the physical problem is far more complex than the idealized problem: for the 
latter explicit formulas for the eigenvalues are known, whereas the eigenvalues of 
the former are given only implicitly as the zeros of an equation involving 4 different 
Bessel functions. 

Example 7.1. Let a = ao be constant and consider the idealized dynamo problem 
(|l.ip . (|1.3|) for which 9 — oo. The eigenvalues of this problem, or equivalently of 
the operator Aoo, can be calculated explicitly (see, e.g., |GK06j ). In fact, since a 
is constant and hence a' = 0, the operator Aoo is given by 

Ax, = ( 'J^ ) , V(Aoo) = 2?(A») © V(Aoo). 

A point A e C is an eigenvalue of Aoo if and only if there is a y — (y±, 2/2)* S T^(Aoo), 
y ^ 0, such that 

(-Ax> - A)yi + a y 2 = 0, 
ao^ooj/i + (-Ax, - X)y 2 = 0. 

Since y% € I?(A 00 ), the first equation yields yi € ©(A 2 ^) and thus the eigenvalue 
equations are equivalent to the relations 

((Aoo + A) 2 - alAoo)yi = 0, a y 2 = {Aoo + X)yi- 

Hence A e o" p (Aoo) if and only if € cTpUAoo + A) 2 — a 2 , Aoo). Now the spectral 
mapping theorem shows that 

&{Aoo) = {-A„(oo) ± a vA^i(oo) : n € N}, 

where A„(oo), n € N, are the eigenvalues of the Bessel operator Aoo with Dirich- 
let boundary conditions (that is, A n (oo) is the n-th zero of A i-> Jj+i/2(vA); see 
Definition l2.il and Remark 13.31 i)). Thus we have 

(7.1) maxa(A 00 ) = -Ai(oo) + \a \ ^1(00) 
and the estimate 

(7.2) maxaf^) = - - M^ + f<^. 

The abstract eigenvalue estimates in Theorem l4.6l and in Theorcm l5.10l combined 
with Remark l5.9l iii). show that, for 9 = 00 and ||a'|| = 0, 

(7.3) maxa(Aco) < aoo — boo = -Ai(oo) + |a | max {|a |, \/Ai(oo) }; 

note that Boo = 5oo and thus 600 = Soo because a is constant (see Proposition ^. II 
and Theorem 15. 3p . 

Comparing (|7.ip and (17.31) , we see that the abstract upper bounds aoo and boo for 
&(<Aoo) in Theorems 14.61 and 15 . 1 01 are sharp for the case ||a|| < ^/Ai(oo). Moreover, 
(|7.2p coincides with the abstract estimate proved in Proposition 15 . 71 for constant a. 
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For physical boundary conditions, even in the case of constant a, only implicit 
formulas are known for the eigenvalues X n (Ai) of Ai (compare |KR80| Section 14, 
in particular, (14.39), (14.41), (14.42)]). 

Example 7.2. Let a = ao be constant and consider the physical dynamo problem 
(|l.ip . (|1.2[) for which 9 = 1. There are no explicit formulas known for the eigenvalues 
of this problem, or equivalently of the operator Ai given by 

Note that here the operators Ai and are given by the same differential expression 
r = — +1(1 + 1) /r 2 , but, unlike the idealized case, they are equipped with different 
boundary conditions (see Definition 13. II) . 

A point A e C is an eigenvalue of Ai if and only if there exists a. y= (2/1,1/2)' G 
L 2 (0, l)8i 2 (0, 1), y ^ 0, with j/i, ^ e ACi oc ((0, 1]), 7% € L 2 (0, 1) for i = 1, 2, such 
that 

(7.4) (-r - A)j/a + a y 2 = 0, 

(7.5) «oTJ/i + (-t - A)t/2 = 0, 
and 

(7-6) y' 1 (l) + ly 1 (l) = 0, 2/2(1) = 0. 

Clearly, the equations (|7.4p . (17.51) are equivalent to the relations 

(7.7) {{t + \) 2 ~alr) yi =0, 

(7.8) a y 2 = (r + A)i/i. 

By Proposition 15. 7[ we already know that for constant a all eigenvalues are con- 
tained in the interval (— oo,ap/4). For A < a§/4, we denote by 

2 I 2 / I 2 \ 2 

ai ai 



the two solutions of the quadratic equation (v + A) 2 — a 2 v = 0. Then we can 
factorize the differential expression in (|7.7|) as 



(7.9) (r + A) 2 - a 2 r = (r - v + (X)) (r - v_ (A)) . 

Since the differential expression r is in limit point case at (see the proof of Pro- 
position [33J , every fundamental system of a differential equation (r — n)x = 
has exactly one solution in £2(0, 1). This and the relation (|7.9[) imply that every 
fundamental system of the differential equation (17. 7|) has exactly two solutions in 
L 2 (0, 1), which may be chosen as the solutions y xl £ L 2 (0, 1) of (r — v±(X))x = 
given by (see (HQ}) 

yt,l(r,X) = J^Vrk±(\)J l+1/2 (rk±(\)) : re [0,1], 



with 



(7.10) fc ± (A) := v^(A) = ^ ± ^ - A. 
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Note that fe+ (A) + k- (A) = ao and fc+(A)fc_(A) = A, that is, fc±(A) are the solutions 
of the quadratic equation k 2 — aok + A = 0. By definition of v±(\) and A;±(A), we 
have (v±(X) + A) 2 = chqV±(X) = alk±(X) 2 . Thus equation (|7.8|) implies that 

yf tl (r,X) = -(r + %±(r,A) = -^(u±(A) + A)j/± (r, A) = fc±(A)y±(r, A). 
ao a 

Therefore a fundamental system y; - } in £2(0, 1) © £2(0, 1) of the system (|7.4I) . 
(|7.5jl is given by 

M) = ( k+ \ x) )vUr,X), y,-(r,A) = (JJ^M)- 

If we take into account the boundary conditions and use relation (|3.10[) . we see 
that there exists a non-zero solution of the boundary eigenvalue problem (|7.4[) . 
(|T31) . (17^)) if and only if 

fc+(A)y+(l) k-(X)y^(l) 

or, equivalently, 

(7.11) Jz_i/ 2 (fc+(A)) J I+1/2 (MA)) - ^+i/ 2 (fc+(A)) Jz_i/ 2 (fc-(A)) - 

with fc±(A) given by (f77TU|) . 

This shows that the eigenvalues of the physical dynamo problem (jl.ip , (|1.2p are 
only given implicitly as the solutions of the equation (|7.11[) . As a consequence, even 
in the simplest case of constant a, it is difficult to obtain any analytic information 
about the eigenvalues of Ai ■ 

The estimates in Theorems 14.61 and 15.101 however, provide the global bound 

maxer(„4;) < min{a;,6;} = bi 

where 

f A 1 (oo) + A 1 (Q / ^M-A^OY , ' ~ " 
b t = I 2 VV 2 ) + Ai(/)|ao| 2 if |aor<Ai(oo), 

[ -Ai(oo) + Kl 2 if |a | 2 >Ai(oo). 

In addition, the estimate in Proposition 15 . 71 for constant a yields the bound 



max(7(*4i) < 




which was already used to establish the eigenvalue relation (|7.11l) . 

Moreover, since Bi = Si for constant a, Proposition 15.131 yields that the number 
ko of positive eigenvalues of Ai coincides with the number of positive eigenvalues 
of the operator — Aoo + a\ or, equivalently, 

# (a(Ai) n (0,oo)) = # (<7(4») n (0,a 2 )); 

if we enumerate the eigenvalues of Ai and of A^ as < Xk Q (Ai) < • • • < Xi(Ai) 
and Ai(oo) < ••• < Afc (oo) < |ao| 2 , respectively, we have the estimate 

< Xj(Ai) < -Aj(oo) + |a | 2 , 3 = l,2,...,fc . 
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Finally, we consider the physical dynamo problem with non-constant helical 
turbulence function a. As shown for the first time in |SG03j . there exist special a- 
profiles which provide dipole-dominated criticality for oscillatory dynamo regimes 
(that is, non-real eigenvalues passing from the left to the right half-plane first for 
dipole modes (I = 1) and later for quadrupole and higher-degree modes (I > 1)). 
Such regimes are of high physical interest due to their close relation to polarity 
reversal processes of the magnetic field (see [SG05J . [SGGX06] ). 

Example 7.3. Consider the physical dynamo problem (jl.ip . (|1.2[) for I = 1 with 
a given by 

(7.12) a(r) =C(-21.46 + 426.41 r 2 - 806.73 r 3 + 392. 28r 4 ), re [0,1], 

where C > is a constant. For C = the eigenvalues are real and they coincide 
with the interlacing sequences of eigenvalues of the Bessel operators — Aoo and —A;. 
According to the numerical computations in |SG03| . if C increases, the largest two 
eigenvalues merge at C — 0.818 and form a complex conjugate pair for 0.818 < 
C < 1.097. In between, at C — 1 this pair crosses the imaginary axis. 

According to the abstract result of Corollary 14. 1 H the largest two eigenvalues 



can only meet if condition (|4.20[) is violated. Indeed, for I = 1, we have 

Ai(0 = vr 2 « 9.87, Ai(oo) « 20.19. 
and, for a given by ([732]) with C = 0.818, 

Hall = max |a(r)| « 17.55, \\a'\\ w max \a'(r)\ w 71.36, 

rG[0,l] re[0,l] 

and hence condition (|4.20j) does not hold since 
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